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BACKGROUND OF THE INVENTION 



The present invention relates to spectroscopic imaging of 



10 heterogeneous light scattering tissue, and more particularly, but not 
exclusively, relates to 7/7 vivo imaging by mapping a fluorescence 
characteristic of the tissue. 

The early detection of disease promises a greater efficacy for 
therapeutic intervention. In recent years, noninvasive techniques have been 
15 developed which have improved the ability to provide a reliable and early 

diagnosis of various afflictions by detecting biochemical changes in the tissue 
of a patient. For example. Magnetic Resonance Imaging (MRI) has 
successfully monitored the relaxation of spin states of paramagnetic nuclei in 
order to provide biomedical imaging and biochemical spectroscopy of tissues. 
20 Unfortunately, the complexity and expense of MRI diagnostics limit its, 
availability — especially as a means of routine monitoring for disease. 

Another powerful analytical technique with an increasing number of 
applications in the biological sciences is fluorescence spectroscopy. 
Applications of fluorescence spectroscopy include biomedical diagnostics, 
25 genetic sequencing, and flow cytometry. As exemplified by U.S. Patent Nos. 
5.421,337 to Richards-Kortum et al. and 5,452,723 to Wu et al.. several 
investigators have suggested various procedures to differentiate diseased 
and normal tissues based on fluorescence emissions through noninvasive 
external measurements or minimally Invasive endoscopic measuring 
30 techniques. Unfortunately, these procedures generally fail to provide a viable 
spatial imaging procedure. One reason imaging based on fluorescence has 
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remained elusive is that meaningful relational measurements of fluorescence 
characteristics from a random, multiply scattering media, such as tissue, are 
difficult to obtain. For example, fluorescent intensity, which is a function of 
the fluorescent compound (or fluorophore) concentration or "uptake" is one 
5 possible candidate for imaging; however, when this property is used in an 
optically dense medium, such as a particulate (cell) suspension, powder, or 
tissue, the local scattering and absorption properties confound measured 
fluorescent intensities. 

Besides intensity, other properties of selected fluorophores such as 
10 fluorescent quantum efficiency and lifetime are also sensitive to the local 
biochemical environment. As used herein, "fluorescent quantum efficiency" 
means the fractional number of fluorescent photons emitted for each 
excitation photon absorbed or the fraction of decay events which result in 
emission of a fluorescent photon. "Fluorescent lifetime," as used herein, is 
15 defined as the mean survival time of the activated fluorophore or the mean 
time between the absorption of an excitation photon and emission of a 
fluorescent photon. Like intensity, measurement of these fluorescence 
characteristics is often limited to well-defined in vitro applications in the 
research laboratory or in flow cytometry where issues such as scattering, 
20 absorption, and changing fluorophore concentrations can be controlled or 
measured. Moreover, these limitations generally preclude meaningful 
fluorescence-based imaging of hidden tissue heterogeneities, such as tumors 
or other diseased tissue regions, which cannot be detected by visual 
inspection. 

25 With the development of techniques to interrogate tissues using 

fluorescence in the near-infrared red (NIR) wavelength regime, noninvasive 
detection of diseased tissues located deep within normal tissues may also be 
possible since NIR excitation and emission light can tfavel significant 
distances to and from the tissue-air interface. U.S. Patent Numbers 

30 5.213,105 to Gratton et al. and 5,353,799 to Chance are cited as further 



background concerning NIR interrogation. As in the case of MRI, x-ray. CT, 
and ultrasound imaging modalities, there is a potential to enhance NIR 
fluorescence imaging techniques with contrast agents. Typically, contrast 
agents for in vivo imaging have depended upon preferential uptake into 
diseased tissue to provide the desired imaging enhancement by absorbing 
the interrogating radiation. The light absorbing tissue provides an enhanced 
spatial variation in measured intensity of the radiation to improve image 
contrast. In the case of a fluorescent contrast agent, the intensity of 
fluorescent light emitted in response to the absorption may provide this 
intensity variation. Generafly, the larger the difference in spatial variation, as 
artificially imposed by a contrast agent, the more improved the reconstructed 
image of interior tissues. Nonetheless, the effectiveness of exogenous 
contrast agents depends greatly upon the selectivity of the agent for the 
tissue region of interest. Unfortunately, targeted and site specific delivery of 
drugs and contrast agents has historically been a limiting factor in both 
therapeutics and diagnostic imaging. Consequently, additional mechanisms 
for inducing contrast that are not dependent solely upon tissue selectively of 
the agent would be,advantageous. 

Thus, a need remains for a technique to noninvasively image multiply 
scattering tissue based on one or more fluorescence characteristics. 
Preferably, this technique includes the implementation of exogenous contrast 
agents with image-enhancing properties beyond preferential absorption of the 
interrogating radiation. The present invention satisfies this need and provides 
other advantages. 
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SUMMARY OF THE INVENTION 

The present invention relates to spectroscopic imaging of 
5 heterogeneous, light scattering nnaterials. Several aspects of the invention 
are novel, nonobvious. and provide various advantages. While the actual 
nature of the invention covered herein can only be determined with reference 
to the claims appended hereto, certain features which are characteristic of the 
present invention are described briefly as follows. 
10 One feature of the present invention is a technique for imaging a 

heterogeneous light scattering material. This technique includes exposing the 
surface of a material to light from a light source and detecting an emission in 
response. A spatial variation of a fluorescence characteristic of the material 
is determined as a function of the emission with a processor. The spatial 
15 variation may be characterized by a set of values representative of the 

fluorescence characteristic as a function of position. An image is generated 
in accordance with the spatial variation that corresponds to the 
heterogeneous composition of the material. This technique may be applied in 
vivo to biologic tissue using external or endoscopic instrumentation to detect 
20 heterogeneities indicative of disease. The technique may include the 
introduction of a fluorescent contrast agent into the material. The 
fluorescence characteristic detected may be fluorescence lifetime, 
fluorescence quantum efRciency, a fluorophore absorption coefficient, 
fluorescent yield (a function of fluorescent quantum efficiency and fluorophore 
25 absorption), or another fluorescence characteristic known to those skilled in 
the art. 

Another feature includes introducing a fluorescent contrast agent into a 
biologic tissue. This contrast agent has a predetermined lifetime and the 
tissue multiply scatters light with a mean time-of-flight between scattering 
30 events. The lifetime and the mean time-of-flight are within a factor of about 
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ten of each other. The tissue is exposed to an excitation light with a 
predetermined time-varying intensity and a light emission is detected from the 
tissue in response to this exposure. An image of the tissue is generated by 
mapping spatial variation of a level of a fluorescence characteristic of the 
5 tissue from the light emission in accordance with a mathematical relationship 
modeling multiple light scattering behavior of the tissue. 

In a further feature, the agent may be selected in accordance with a 
predetermined relationship between degree of image contrast and at least 
one of fluorescence yield or the fluorescence lifetime. Preferably, the lifetime 
10 is in a range of about 0.1 to 10 nanoseconds (ns). A more preferred range is 
0.5 to 5 ns. A still more preferred range is about 0.2 to 2 ns. A most 
preferred value for the lifetime is about 1 ns. 

An additional feature includes evaluating ability of a number of 
fluorescent agents to provide image contrast between different tissue types. 
15 This evaluation includes determining a relationship between degree of image 
contrast and at least one of fluorescence lifetime or fluorescence yield of the 
agent. One of the agents is selected based on the evaluation. The selected 
agent is provided for introduction into a biologic tissue to enhance imaging 
performed in accordance with a mathematical expression modeling the 
20 behavior of multiply scattered light traveling through the tissue. 

In still another feature, a biologic tissue is exposed to a first excitation 
light and a first emission is detected from the tissue in response to the first 
excitation light. A fluorescent contrast agent is introduced into the tissue after 
this detection and the tissue is then exposed to a second excitation light. A 
25 second emission is sensed in response to the second excitation light. Data 
corresponding to the first emission is compared with data corresponding to 
the second emission to evaluate contrast provided by the agent. Contrast 
may be determined as a function of at least one of fluorescence lifetime, 
fluorescence yield, or quantum efficiency. For a frequency domain form of 
30 this evaluation, the image contrast may be evaluated in terms of phase 
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contrast or modulation contrast Moreover, the wavelength of the first 
excitation light may be selected to be generally the same as the wavelength 
of the fluorescent light emitted by the agent in response to the second 
excitation light. 

5 Accordingly, it is one object of the present invention to map a 

fluorescent property of a light scattering material that varies with the 
heterogeneous composition of the material to generate a corresponding 
image. 

It is another object of the present invention to provide a spectroscopic 
10 technique for noninvasively monitoring fluorescent properties of hidden tissue 
volumes in a living organism and to monitor selected metabolites of an 
organism in vivo. 

Yet another object is to provide a technique to select and design 
fluorescent contrast agents who improve contrast for photon migration based 
15 imaging. This technique may include the selection of contrast enhancing 
properties that are not solely dependent upon uptake. 

Further objects, features, aspects, benefits, and advantages of the 
present invention will become apparent from the drawings and description 
contained herein. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a schematic illustration of a system of one embodiment of the 
present invention. 

5 FIG. 2 is a flow chart of a process utilizing the system of FIG. 1, 

FIG. 3 is a schematic representation of a tissue phantom arrangement 
used to demonstrate various aspects of the present invention. 

FIGS. 4-7 graphically depict selected properties of equations used in 
the present invention. 
10 FIGS. 8 and 9 graphically depict convergence of simulated 

determinations of the spatial variation of fluorescent yield and lifetime, 
respectively, utilizing one embodiment of the present invention. 

FIGS. 10-14 are computer-generated gray scale images obtained from 
experimental examples 1-3 of the present invention. 
15 FIG. 15 is a schematic illustration of a system of an alternative 

embodiment of the present invention. 

Fig. 16 is a schematic view of another tissue phantom arrangement of 
the present invention. 

Fig. 17A is a graph of simulated measurements of phase-shift (vertical 
20 axis) as a function of angular detector position (horizontal axis) to compare a 
contrast agent with a 1 nanosecond (ns) lifetime (solid line with open 
symbols) to a contrast agent with a 1 millisecond (ms) lifetime (dashed line 
with closed symbols); where the contrast agents are selective to a 
heterogeneity and the different line styles each correspond 10 millimeter 
25 (mm), 20 mm. 30 mm, 40 mm, and 50 mm positions of the heterogeneity. 
Fig. 17B is a graph of simulated measurements of modulation 
amplitude (vertical axis) as a function of angular detector position (horizontal 
axis) to compare a contrast agent with a 1 ns lifetime (solid line with open 
symbols) to a contrast agent with a 1 ms lifetime (dashed line with closed 
30 symbols); where the contrast agents are selective to a heterogeneity and the 
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different line styles each con-espond 10 mm. 20 mm. 30 mm, 40 mm, and 50 
mm positions of the heterogeneity. 

Fig. 18A is a graph of simulated measurements of phase contrast 
(vertical axis) as a function of detector position (horizontal axis) and location 
5 of a heterogeneity (different line styles corresponding to 10 mm, 20 mm. 
30mm, 40mm, and 50mm positions); where the heterogeneity contains a 
light-emitting contrast agent having an uptake into the heterogeneity of 100:1 
and a lifetime of 1 ns. 

Fig. 18B is a graph of simulated measurements of phase contrast 
10 (vertical axis) as a function of detector position (horizontal axis) and location 
of a heterogeneity (different line styles corresponding to 10 mm, 20 mm, 
30mm, 40mm, and 50mm positions); where the heterogeneity contains a 
light-emitting contrast agent having an uptake into the heterogeneity of 100:1 
and a lifetime of 1 ms. 
15 Fig. 18C is a graph of simulated measurements of modulation contrast 

(vertical axis) as a function of detector position (horizontal axis) and location 
of a heterogeneity (different line styles corresponding to 10 mm. 20 mm, 
30mm, 40mm,. and 50mm positions); where the heterogeneity contains a 
light-emitting contrast agent having an uptake into the heterogeneity of 100:1 
20 and a lifetime of 1 ns. 

Fig. 18D is a graph of simulated measurements of modulation contrast 
(vertical axis) as a function of detector position (horizontal axis) and location 
of a heterogeneity (different line styles corresponding to 10 mm, 20 mm, 
30mm. 40mm, and 50mm positions); where the heterogeneity contains a 
25 light-emitting contrast agent having an uptake into the heterogeneity of 100:1 
and a lifetime of 1 ms. 

Fig. 19A is a graph of experimental measurements of phase contrast 
(vertical axis) of emission light as a function of detector position (horizontal 
axis) and heterogeneity location (different line styles corresponding to 10 mm. 
30 20 mm, 30 mm, and 40 mm) for a 100: 1 uptake of an ICG contrast agent into 
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the heterogeneity. 

Fig. 19B is a graph of experimental measurements of phase contrast 
(vertical axis) of emission light as a function of detector position (horizontal 
axis) and heterogeneity location (different line styles corresponding to 10 mm. 
5 20 mm. 30 mm. and 40 mm) for a 100:1 uptake of an Ru(bpy)3^* contrast 
agent into the heterogeneity. 

Fig. 20A is a graph comparing absorption measurements (open 
symbols) and fluorescent measurements (closed symbols) in terms of phase 
contrast (vertical axis) versus heterogeneity (object)--position in centimeters 
10 (cm) (horizontal axis) for an ICG contrast agent at modulation frequencies of 
80 and 160 megahertz (MHz) (different symbol shapes). 

Fig. 20B is a graph comparing absorption measurements (open 
symbols) and fluorescent measurements (closed symbols) in terms of 
modulation contrast (vertical axis) versus heterogeneity position (horizontal 
15 axis) for an ICG contrast agent at modulation frequencies of 80 and 160 
megahertz (MHz) (different symbol shapes). 

Fig. 20C is a graph comparing absorption measurements (open 
symbols) and fluorescent measurements (closed symbols) in terms of phase 
contrast (vertical axis) versus heterogeneity position (horizontal axis) for an 
20 DTTCI contrast agent at modulation frequencies of 80 and 160 megahertz 
(MHz) (different symbol shapes). 

Fig. 20D is a graph comparing absorption measurements (open 
symbols) and fluorescent measurements (closed symbols) in terms of 
modulation contrast (vertical axis) versus heterogeneity position (horizontal 
25 axis) for an DTTCI contrast agent at modulation frequencies of 80 and 160 
megahertz (MHz) (different symbol shapes). 

Figs. 21A-21D are computer-generated gray scale images formed from 
experimental measurements depicting spatial variation for a tissue phantom 
including separated ICG and DTTCI heterogeneities in terms of modulation 
30 phase-shift, AC modulation amplitude, average DC intensity, and modulation 
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ratio (AC/DC), respectively. 

Figs. 22A-22D are computer-generated gray scale images depicting 
spatial variation for in vivo tissue Imaging of a dog treated with an ICG 
contrast agent in temns of modulation phase-shift, modulation ratio (AC/DC), 
5 average DC intensity, and AC modulation amplitude, respectively. 
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DESCRIPTION OF PREFERRED EMBODIMENTS 

For the purposes of promoting an understanding of the principles of the 
invention, reference will now be nriade to the embodiments illustrated in the 
5 drawings and specific language will be used to describe the same. It will 
nevertheless be understood that no limitation of the scope of the invention is 
thereby intended. Any alterations and further modifications in the described 
techniques, methods, systems, and devices; and any further applications of the 
principles of the invention as described herein are contemplated as would normally 
10 occur to one skilled in the art to which the invention relates. 

FIG. 1 depicts system 110 of the present invention for fluorescent 
imaging of tissue 100. Tissue 100 has surface 101 and a heterogeneous 
composition as represented by regions 102, 103 underlying surface 101. 
Heterogeneities 102. 103 are generally not detectable by visual inspection of 
15 surface 101. 

System 110 includes modulated light source 120 to supply an intensity- 
modulated excitation light of predetermined frequency and wavelength to 
tissue 100 via optic fiber 123. Preferably, source 120 is a laser diode of 
conventional design with a modulated output in the 1-500 MHz frequency 
20 range and a monochromatic output in the 100 to 1000 nanometer (nm) 

wavelength range. The specific wavelength is selected to excite a designated 
fluorophore in tissue 100, Beam splitter 126 may be employed to direct a 
small portion of the excitation signal to reference sensor 128 for processing 
purposes. 

25 System 110 also includes detection subsystem 140 which has optic 

fibers 143 to detect photons emitted from tissue 100 from a number of 
corresponding detection sites. Subsystem 140 includes one or more 
emission sensors 148. Detection subsystem 140 also includes an 
interference filter to obtain a selected emission wavelength corresponding to 

30 emission of a designated fluorophore in tissue 100. In one embodiment. 
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subsystem 140 includes a single sensor 148 and the signals from fibers 143 
are multiplexed. Preferably, sensors 128, 148 are Photo-multiplier Tubes 
(PMTs) or photodiode detectors but other sensor varieties, such as image 
intensifiers and charge-coupled devices, are also contemplated. 

Sensors 128, 148 and source 120 are operatively coupled to 
heterodyne subsystem 130. Subsystem 130 is configured to obtain 
infomnation about the phase, AC, and DC intensity of light detected with 
sensor 128 relative to light detected with the sensor 148 using conventional 
heterodyning techniques. In one embodiment, heterodyne subsystem 130 
includes a signal synthesizer phase-locked to the repetition rate of a laser 
used for source 120. For this embodiment, subsystem 130 includes an 
amplifier to gain modulate sensors 128, 148 at a harmonic of a laser 



repetition rate (when a pulsed laser is used) or at the modulation frequency 
(when a modulated laser diode is used) plus an offset to provide the desired 
heterodyning. In one variation of this embodiment, an 80 MHz pulsed laser 
repetition rate is divided down to 1 0 MHz and input to the synthesizer, and a 
heterodyning offset of 100 kHz is input to the amplifiers for sensors 128, 148. 

Sensors 128, 148 are operatively coupled to processor 160. 
Processor 160 includes input/control device 162, output device 164, and 
memory 166. Processor 160 may be an electronic circuit comprised of one or 
more components. Similarly, processor 160 may be comprised of digital 
circuitry, analog circuitry, or both. Also, processor 160 may be 
programmable, an integrated state machine, or a hybrid combination thereof. 
Preferably, input device 162 is a keyboard or input control of a conventional 
variety, and output device 166 is a Cathode Ray Tube (CRT) based video 
display, printer, or other image display system known to those skilled in the 
art. Memory 166 is preferably of the electronic (e.g. solid state), magnetic, 
or optical variety of the type readily available for use with electronic 
controllers or processors. Furthermore. Memory 166 may include an optical 
disk memory (CD), electromagnetic hard or floppy disk media, or a 
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combination of these, 

FIG. 2 depicts one nnode of operation of system 1 10 as process 210. 
Process 210 includes mapping the spatial variation of fluorescence yield and 
lifetime with processor 160 and generating an image signal in accordance 
5 with the map. Output device 164 is configured to display an image in 
response to the image signal. Process 210 begins by introducing a 
fluorescent contrast agent into tissue 100 in stage 212. This agent provides a 
source of fluorescent emission for detection by subsystem 240. The 
configuration of the modulated light source 120. heterodyne subsystem 130, 
10 and detection subsystem 140 is designed to accommodate the excitation and 
emission properties of the selected fluorescent agent. In other embodiments, 
endogenous fluorophores may be alternatively or additionally employed and 
system 110 adjusted accordingly. 

In stage 214. light source 120 configured according to the selected 
15 fluorophore excites tissue 100, In stage 216. the phase. Oobs- and log of AC 
intensity, Mobs, of the emission at each detection site "i" relative to the 
excitation light from source 120 are determined at the heterodyne (or offset) 
frequency. For "Di" number of detection sites, the detected or observed 
phase and AC intensity are indexed by "i" using the following notation: 
20 (9obs)i and (Mobs)i. respectively. Processor 160 stores the relative phase 
and AC intensity information in memory 166. 

In stage 218, a two dimensional grid is established for an area of 
tissue 100 selected for imaging, and a matrix of grid points is established and 
indexed by "j", A uniform seed value for the fluorescent yield, yj = (r|Max->m^j' 
25 and the fluorescent lifetime. (x)j. at each grid point j is assigned. These 
values are an initial homogeneous guess of the yield and lifetime values, 
which are modified in later stages. The term "ri" is the quantum efficiency of 
the fluorophore which varies with the environment of the surrounding of the 
fluorophore. The term "Max— ►m * absorption coefficient for the 
30 fluorophore and is the product of the extinction coefficient of the fluorophore 
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based on the natural log and the concentration of the fluorophore. As a 
result, the yield, y =nMax-^m' influenced by the sun-ounding metabolism 
and the uptake of the fluorophore. The uptake of certain known fluorophores 
vary with the type and condition of host tissue, providing another fluorescence 
characteristic useful to detect disease. The contrast provided by these 
properties is largely independent of fluorophore concentration. The initial 
estimate of fluorescent yield and lifetime are stored in memory 166 by 
processor 160 for later use. 

After establishing this initial estimate of the fluorescence characteristics 
of yield. HMax-^m* lifetime, t. processing loop 220 is entered in stage 
230. Preferably, the stages of processing loop 220 are executed by 
processor 160 via preprogrammed software, dedicated hardware, or a 
combination of both as appropriate. To aid in understanding various 
mathematical aspects of process 210 and loop 220, the following table of 
selected variables is listed: 

c velocity of light in a vacuum; 

D(r) optical diffusion coefficient; 

Oi number of detection sites; 

f modulation frequency; 

I identity matrix; 

i detection site index; 

J Jacobian matrix relating the sensitivity at each grid point, 



J, 



j, to the response at each detection site; 
grid point index; 

individual elements of the Jacobian matrix J; 



k 



source index; 

log of AC intensity of modulated fluorescent light position; 
index to multiple modulation frequencies; 
average index of refraction; 



M 



m 



n 
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r 



position (in two or three dimensions); 
Sk number of modulated light sources; 

S(r. w) source term for the modulated light at position r and 

frequency cu; 



Greek 

merit function representing the least squares error; 
0x(''» ^) complex number representing photon flux in the 

Frequency domain at position r and frequency co; 
10 n quantum efficiency of fluorescent probe or dye; 

Ma average absorption coefficient; 

Map^ absorption coefRcient of the fluorescence light by both 

the nonfluorescing chromophores and fluorophore; 
Max absorption coefficient of the excitation light by both the 

15 nonfluorescing chromophores and fluorophore; 

Max— >c adsorption coefficient due to nonfluorescing 
chromophores; 

Max— >m adsorption coefficient of excitation light by fluorophores; 
M's effective scatting coefficient; 

20 9 phase-shift of one modulated light wave to another; 

X lifetime of activated probe or dye at location r; 

uj angular modulation frequency, given by 2nt 



Subscripts 

25 obs observed or experimental data; 

X excitation light; and 

m fluorescence or emission light. 



In stage 230, phase and relative AC intensity at each detection site "i" 
30 is calculated as a function of the initial estimates of yield and lifetime for each 
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grid point j. The caiculated phase and intensity are represented at each 
detection site i as (em)i and (Mm)j, respectively. The values for (Gm)} and 
(Mm)i are determined using the diffusion equation approximation of the 
radiative transport equation. The diffusion equation approximation describes 
5 the spatial and temporal transport of light in tissues or multiply scattering 
media. A coupled frequency domain diffusion equation can be used to 
predict the excitation and emission fluence rates. (Dx(r, cu) and <t>m (r. c*)). 
respectively, at any location r within the selected grid of tissue 1 00 via 
equations (1) and (2): 

10 

V«[Dx(r)V<Dx(r.(o)HMax('')+iw/Cn]<l>x(r.u)) +Sx(r,ca) = 0 (1) 
V»[Dm(r)V<I>m(r.cj)J-[Mam('')+'^^Cn]<I>m('".w) +Sm(r.a)) = 0 (2) 

15 

The source term for the excitation light Sx(r.CL)) is due to the sinusoidally 
modulated light at an angular frequency w = 2Trf where f is typically in the 
MHz frequency range. The first term in both of the dHfusion equations (1) and 
(2) represents the diffusive or "random-walk" transport of multiply scattered 
20 light where Dx,rn 'S the optical diffusion coefficient of equation (3) as follows: 

Dx.m=l3(Pax.m*M'sx.m)l"'' (3) 

and Pa ^nd m ^ are the absorption and isotropic scattering coefficients, 
25 respectively, for tissue 100, the medium of interest. Multiple scattering of light 
occurs when p g » pg; where pg indicates the ability to absorb light and p'g 
indicates the ability to scatter light for a given material at a given wavelength. 
As used herein, "multiply scattered light" r^rers to light that travels at least five 
(5) times the mean isotropic scattering length, defined as l/pg. 
30 Because these optical properties are dependent on the wavelength of 
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light, the coefficients generally differ for the excitation light from source 120 
(subscript x) and fluorescent emission detected with subsystem 140 
(subscript m). The total absorption coefficient at the excitation wave*ength, 
jjg^. is due to contributions from nonfluorescing chromophores as well as 
5 from fluorophores responsive to the excitation wavelength. The total 

absorption coefficient is given by the sum of absorption coefficients due to 
nonfluorescing chromophores, Max— *^c* fluorophores Pax-H^rrr G©n®rally 
it may be assumed that the absorption experienced at the fluorescent 
wavelength is due primarily to nonfluorescing chromophores. The velocity of 
10 light in tissue is Cr,=c/n where n is the average index of refraction. The source 
term for the fluorescent emission is dependent on the excitation light fluence, 
<Dx(r. CO) and is given by equation (4) as follows: 

Sm(r,a))=nMax_,rn ('")^x(^ co)[(1-iojx(r))/(1+a)2x(r)2)] (4) 

15 

This term arises from the Fourier transform of the fluorescence decay term in 
the time domain following an incident pulse of excitation light where: x Is the 
fluorophore lifetime, q is the quantum efficiency, and the absorption 
coefficient. Max— *^m* product of the extinction coefficient based on 

20 natural log and the concentration of the fluorophore in the ground state. As 
previously indicated, the combined product. HMax— >m' '® termed the 
fluorescent yield, y, ^ind Is proportional to the generated fluorescence fluence. 
Substitution of equation (4) into equation (2) facilitates determination of 
for each grid point "j." The solution of the diffusion equations (1) and (2) for 

25 the two-dimensional area defined by the grid points "j" may be readily 

extended to three dimensions to estimate spatial variation of one or more 
fluorescence characteristics in a selected volume with "r" corresponding to 
position in three dimensions. 

Both diffusion equations (1) and (2) are linear complex elliptic 

30 equations that can be solved as boundary value problems for the complex 
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quantities O^Cr, co) and Or^Cr co). This solution employs the nnethod of finite 
differences to create corresponding finite difference equations. These 
difference equations are utilized to obtain an approximate solution at each 
grid point, j. This method of solution is described in other contexts in Fulton 
5 et al.. Multigrid Method for Elliptic Problems, A Review , 1 14 American 

Meteorological Society pp. 943-59 (May 1986); and B.W. Pogue et al.. Initial 
Assessment of a Simple System for Frequency Domain Diffuse Optical 
Tomography , 40 Physics in Medicine and Biology pp. 1709-1 729 (1995). One 
preferred method of performing this solution is with the MUDPACK routines 

10 described in Adams, J.C., MUDPACK: Multigrid Portable Fortran Software for 
the Efficient Solution of Linear Elliptic Partial Differential Equations , 34 App. 
Math Comp. p. 133 (1989). For the solution of the diffusion equations, it is 
assumed that Opn,x(^^) = 0 on the surface 101 of tissue 100 which is known 
as the zero fluence boundary condition. It should be recognized that other 

15 boundary conditions may be selected and the method of solution varied 
accordingly. 

The diffusion equations (1) and (2) may be solved for a complex 
number for Om at each grid point, j. The detected signal at the surface is 
proportional to the normal component of the gradient of the photon fluence. 

20 To approximate the signal at detector site "i" located on surface 101 of tissue 
100, the cDppi value at an internal grid point closest to the site is selected which 
follows from the relationship that the normal component of the photon fluence 
gradient is proportional to <t>^ just inside surface 101 . The calculated phase- 
lag, 9m, and the log of AC intensity. Mm. at the detection sites "Di" are 

25 calculated from the imaginary and real parts of the complex <Dm value with 
respect to the phase and the AC intensity of source 120. 

The diffusion equations (1) and (2) provide insight into the sensitivity of 
changing the fluorescent optical properties of tissue 100 on 9^ and Mm 
measured at the detector sites 1. This insight results from a series of 

30 calculations while fixing various parameters of the diffusion equations (1) and 
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(2). These calculations assume circular tissue phantom 300 with an 
embedded, heterogeneity 302 hidden in phantom background 303 as 
illustrated In FIG. 3. A two-dimensional grid is established for phantom 300 
and may easily be expanded to three dimensions. Under these simulated 
5 conditions, a large value is assigned to absorption coefficients for both 
excitation and fluorescent light at all grid points outside the simulated tissue 
phantom. The four sources S1-S4 of FIG. 3 (Sk = 4) are simulated by 
assigning an arbitrary complex number at a grid point near the surface closest 
to each source. The twenty detection sites D1-D20 of FIG. 3 (Di = 20) are 
10 simulated by using the calculated values detemiined from cd^ at the grid 
point "j" closest to the detection site. The simulated solutions to diffusion 
equations (1) and (2) were obtained in two dimensions for a 65 x 65 grid 
covering a 1 00 mm diameter circular tissue phantom 300 with a circular, 
embedded heterogeneity having a 30mm diameter and located at the center 
15 of the tissue phantom 300 (this location differs slightly from the configuration 
of heterogeneity 302 of FIG. 3). The simulated measurements of fluorescent 
phase-shift and AC intensity are reported for 20, equally spaced, 
circumferentially located detection sites D1-D20. The modulation frequency, 
f, was set equal to 150 MHz. The optical properties of the heterogeneity and 
20 the background are shown in Table ^ as follows: 
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T 


frequency 










background 


background 


(MHz) 






(mm-^) 






(ns) 






0.0 


1.0 


0.0 


1.0 X 10-5 


1.0 


150.0 
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In order to evaluate the influence of nMax-*m' ®m Mm were 




20 



computed at each detection site D1-D20 as the value of nMax-*m 
heterogeneity increased from 10"^ mm"'' to 10"^ mm"^ and as HMax— ►m 
the background 303 was maintained constant. The lifetime, x, was set equal 
to 1 ns for both the object and the background causing contrast due to 
5 differences in n|Jax->m- P'^*^ ^m and M^^ are shown in FIGS. 4 and 5 
respectively for one active source S1 . As nMax-*m heterogeneity 102 
increases to higher values, the AC intensity approaches an upper limit similar 
to what is expected in dilute non-scattering solutions. FIG. 5 shows how the 
fluorescent phase-shift, 9^^, decreases as the absorption coefficient due to 
10 the fluorophore, Max— ..m decreased 10 to 100 times the background. From 
these simulations, Mp^ appears to be directly dependent upon changes in 
HMax-^m simulated tissue heterogeneity 102 whereas 9^ is indirectly 
dependent on HMax— ►m changes in photon migration. 

In order to evaluate the influence of t, 9^. and Mm were calculated at 
15 each detection site D1-D20 as the values of t in the heterogeneity varied from 
10"'' ns to 103 ns and the value of x in the background was held at 1 ns. The 
background nMax-*m lO'^mm-^ and HMax-^m ^^''the 

heterogeneity was set to lO-^mm*^. As shown in Fig. 6, the detected AC 
intensity increases as x decreases. Fig. 7 illustrates the values of the 
20 fluorescent phase-shift at each detection site as the lifetime of the 

heterogeneity is changed from 0.1 ns to 1000 ns. At a given modulation 
frequency (150 MHz in this calculation). 9m first increases, reaches a 
maximum and then subsequently decreases as x is increased from 0.1 ns to 
1000 ns. Therefore, both 9m and Mm at each detection site D1-D20 appear to 
25 be directly influenced by the value of lifetime in the heterogeneity. 

Referring back to FIG. 2. in stage 240. the calculated emission phase 
and intensity, (9m)i and (Mm)i. are compared to the measured emission 
phase and intensity, (9obs)i and (Mobs)i' breach detection site "i" to identify 
a difference or "error" between the measured and calculated values. 
30 Because (hMax— #m^j *^Pacts (Mm)j . this comparison is posed in the form of 
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the merit function of equation (5) as follows: 

Sk Di 

Xp2 = (1/Sk) Z (1/Di) 2 I((Mobs)i -(Mm)iVc^Mj2 (5) 
5 k=1 i=1 

where om is the typical standard deviation of noise In Mm, taken to be 0.01 ; 
Sk = number of excitation source sites indexed by k; and Di = number of 
detection sites indexed by i. The goal of the algorithm is to minimize by 
10 appropriate updates of (nMax—m^j- ^^^^ ^n initial update of (nPax— m^j 

another merit function in terms of (T)j participates in the comparison of stage 
240. This merit function, jg presented as equation (6) as follows: 

Sk Di 

15 Xt^ = (1/Sk) Z (1/Di) Z [((Mobs)i-(Mm)i)/ ou ]2+I((eobs)i-(em)i)/cre]2 (6) 
k=1 i=1 

where oq is the typical standard deviation of noise in (9^)1. taken to be 1 
degree; Sk = number of excitation source sites indexed to k; and Di = 

20 number of detection sites indexed to i. Since the lifetime influences both 
(6m)i ^^'^ (Mm)i. the phase and AC intensity are used in equation (6), 

After the comparison of stage 240 is performed by calculating the merit 
functions X(j2 and Xt^ control flows to conditional 250 to test whether the 
comparison of the observed values, (Oobs)! and (Mobs)i. to the calculated 

25 values (9ni)i and (Mfn)j via the merit functions meets a selected convergence 
criteria. This criteria corresponds to the degree of tolerable error in 
determining the yield and lifetime values. For one embodiment, convergence 
is achieved when any of the following three quantities, (i) x^, (ii) change in x^ 
in successive iterations of loop 220, or (iii) relative change in x^ in successive 

30 iterations of loop 220 is less than a predetermined threshold value of 
1.0 X In other embodiments a different comparison calculation and 
associated conditional may be employed as would occur to one skilled in the 
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art. If conditional 250 is satisfied, control flows to stage 270 and loop 220 is 
exited; however, if the criteria is not satisfied, execution of loop 220 continues 
in stage 260. 

In stage 260. the yield. (y)j= (nMax-rn)]' lifetinne. (T)j. for each grid 
5 point j is updated so that these values may reach the nninimum error 

corresponding to the comparison stage 240 and conditional 250 test. In order 
to update these values. Jacobian matrices are used which describe the 
sensitivity of the response at each detection position i to changes in (y)j= 
(HMax-^m^j- lifetime. (x)j, at each grid point, j. Three Jacobian matrices 
10 are employed: 

J(M.nMax-.m): ^(M.x): JO.t); 

where, the elements Ji.j of these Jacobian matrices are given by Ji.j = 
15 [5Mi/(5(niJax-.m)jl: ^'-i = [5Mj/aTj]; and Ji.j = [dQ/dxjl respectively. These 

elements may be calculated by solving the diffusion (1) and (2) four times for 
each grid point, j to obtain j and e^.i calculated with (t)] and (T+aT)j and 

with (HMax— m)j (HPax-^m "*"*^nMax-..m^j- ^^^^ squares 
minimization, the update to yield and lifetime is calculated. In one preferred 

20 embodiment, this updating algorithm is adapted from an algorithm used to 
reconstruct images obtained by electrical impedance tomography like the 
algorithm suggested by Yorkey, et aL , Comparing reconstruction Algorithms 
for Electrical Impedance Tomography . 34 Transactions in Biomedical 
Engineering pp. 843-52 (1987). The Jacobian Matrices are used to solve for 

25 update vectors, [AnMax-+m ^ estimated yield and lifetime vectors. 

[TMax— m ^ respectively. These vectors are of a dimension 

corresponding to the number of grid points. At each iteration through loop 
220, the following Jacobian equations (7) and (8) are solved to determine the 
update for the estimated yield and lifetime vectors: 
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$ ^ (Mm^ - Mm) 



[ 
[ 



j(M.T)'^j(M.T) j(e.T)''' j(e. T) 



5 



(7) 



(8) 



Mfnob, '^m 3re the observed and calculated vectors of the log of AC 
intensity at each of the i detection sites, respectively, '^m^^ ^re the 

observed and calculated vectors of the phase lag at each of the i detection 
15 sites, respectively. Due to the ill-conditioned nature of the Jacobian matrices, 
the terms A-jl or A2I are added as part of a Marquardt minimization scheme 
where I is an identity matrix. The parameters A-j or A2 are adjusted via a 
Maquardt-Levenberg type algorithm of the type disclosed in Press et al., 
Numerical Recipes: The Art of Scientific Computing , (Cambridge University 

20 Press, 1992). Conventional numerical methods are employed to solve the 
simultaneous linear algebraic equations resulting from the Jacobian matrix 
equations (7) and (8). The Jacobian matrices are re-caiculated at each 
iteration through loop 220. It has been found that equations (7) and (8) 
provide a way to select appropriate changes to the yield and lifetime 

25 estimates; however, other numerical approaches to recursively iterate to 
acceptable estimates as would occur to one skilled in the art are also 
contemplated. Once the update is complete, control returns to stage 230. 

If the convergence criteria is satisfied in conditional 250. then 
estimation of yield and lifetime for the grid points has reached an acceptable 

30 minimum and control flows to stage 270. In stage 270 an image signal Is 
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generated by processor 160 from the spatial variation of the yield and/or 
lifetime fluorescence characteristics. This image signal is sent to output 
device 164, which displays an image in response. Because the fluorescence 
characteristics of yield and lifetime typically vary with the biologic environment 
5 of the fluorophore. this image is generally indicative of tissue variation and 
offers the capability to detect heterogeneities 102. 103. For example, laser 
diodes capable of supplying Near infrared (NIR) light that can penetrate tissue 
several centimeters, and fluorescent contrast agents responsive to NIR light 
may be used to provide a viable imaging system. In one embodiment, this 
10 system is adapted for use with an endoscope. 

Besides yield and lifetime, the spatial variation of other fluorescence 
characteristics useful to distinguish diseased tissues may be mapped using 
the diffusion equations (1) and (2). Such alternative fluorescence 
characteristics include, but are not limited to. quantum efficiency ri and/or 
15 fluorescent absorption coefficient Max— >m determined as separate properties 
independent of the yield product. 

It should be appreciated that imaging in accordance with the present 
invention, such as process 210. includes exposing biologic tissue at the 
tissue-air interface to an excitation light and detecting the light which has 
20 propagated to a detector located some distance away from the source on the 
air-tissue interface. The time-dependent propagation characteristics of 
multiply scattered light emitted in response to this exposure are measured. 
As described in connection with process 210, an intensity-modulated light 
source may be employed for frequency-domain measurements. The 
25 propagating wave of intensity-modulated light is amplitude attenuated and 
phase-shifted relative to the excitation light owing to the spatial distribution of 
fluorescence properties. From exterior measurements of phase-delay and 
amplitude modulation, interior fluorescence properties are determined using a 
mathematical relationship that models the multiple light scattering behavior of 
30 the tissue, such as the diffusion equations (1) and (2). These fluorescence 




25 



properties may be mapped to provide a corresponding interior image of the 
tissue, facilitating the identification of hidden heterogeneities. 

In an alternative embodiment, measurements may be made in the time 
domain. For this embodiment, a pulse of light may be launched at the air- 
5 tissue interface, which is broadened during its propagation in tissues due to 
the spatial variation of fluorescence properties within the tissue. The 
broadened pulse emitted from the air-tissue interface is measured. For this 
embodiment, the diffusion equation in the time-domain form, or such other 
mathematical relationship characterizing multiply scattered light propagation 
10 through the tissue may be utilized to calculate the ffuorescence 
characteristics as would occur to those skilled in the art. These 
characteristics may then be mapped to generate a corresponding image in 
the manner described in connection with process 210. 

Both the frequency and time domain approaches account for the time 
15 propagation of light through the tissue due to multiple scattering events. For 
a given photon, the travel time through a multiple scattering media increases 
with the number of collisions or "scattering events", which corresponds toB 
longer scattering path. This travel time is known as the "time-of-flighf. 
Typically, time-of-flight is on the order of a fraction of a nanosecond to a few 
20 nanoseconds in biologic tissue. For the usual case of many photons each 
traveling along different scattering paths, a mean "time-of-flight" of the 
photons may be determined from the frequency or time domain 
measurements. These time-based measurements are utilized with the 
corresponding mathematical model to map the fluorescence characteristics. 
25 The fluorescence characteristic map provides an image of tissue that 

may be based only on intrinsic fluorophores in the tissue or enhanced by 
introduction of a contrast agent that is selective to tissue volumes of interest. 
This contrast agent may absorb radiation as in the case of contrast agents for 
x-ray and CT imaging to provide a corresponding darkening of the image 
30 regions for the tissue volumes of interest. Unfortunately, the contrast 
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provided through selective absorption is limited. Accordingly, in another 
embodiment of the present invention, a technique to select exogenous 
contrast agents which augment the conventional contrast mechanisms is 
provided. It has been discovered that fluorescence properties that change 
5 with the local biochemical environment often provide greater contrast for 

reconstruction of diseased tissue volumes that can be afforded by absorption- 
based contrast alone. Among the properties that offer this local environment 
contrast mechanism, are fluorescence lifetime t, i.e., the mean time between 
the absorption of an excitation photon and the emission of a fluorescent 
10 photon; fluorescence quanturh efficiency n. i.e.. the number of fluorescent 
photons emitted per excitation photon absorbed; and fluorescence quantum 
yield y. 

It has been discovered that fluorophore contrast agents having a 
fluorescence lifetime within an order of magnitude - or factor of ten (10) — of 

15 photon "time-of-flights" of the tissue being interrogated are surprisingly 

advantageous in providing contrast for photon migration imaging. One way of 
utilizing this surprising advantage is to select an agent with a fluorescence 
lifetime within a factor of ten (10) of the mean time-of-flight predicted for the 
tissue to be imaged. Typically, by applying this principle, a preferred range 

20 for the contrast agent lifetime of about 0.1 to 10 ns results. More preferably, 
the range for the fluorescence lifetime of the contrast agent is within a range 
of about 0.5 to 5 ns. A still more preferred range for fluorescence lifetime of 
the agent is about 0.2 to about 2 ns. A most preferred value for the lifetime is 
about 1 ns. 

25 It has also been discovered that fluorescence characteristics may 

influence the resolution of measurements of the detected light emission. For 
example, in the frequency domain, it has been found that the amplitude of the 
intensity-modulated fluorescent light emanating from a hidden heterogeneity 
containing the agent generally increases with quantum yield y or quantum 

30 efficiency r\. Further, as fluorescence lifetime t within the heterogeneity 
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increases relative to its surroundings, the phase contrast increases. 
Conversely, the amplitude of the detected intensity-modulated light decreases 
with increasing fluorescence lifetime t within the heterogeneity relative to its 
surroundings. Through these discoveries, a fluorescent agent may be 
5 selected or formulated to provide a desired measurement resolution and 
fluorescence lifetime contrast suitable for photon migration interrogation of a 
heterogeneous arrangement of tissue. These discoveries are further 
described in connection with Examples 4-7 at the end of this description. 

Generally, as Examples 4-7 illustrate, the selection or formulation of a 
10 suitable contrast agent is performed by determining the relationship between 
image contrast and fluorescent properties such as lifetime, yield, or quantum 
efficiency a function of the location of a heterogeneity selective to the given 
contrast agent. These relationships may be evaluated for a number of 
different agents to select a preferred agent for a given contrast problem. For 
15 frequency domain based evaluation, image contrast may be characterized in 
terms of phase shift variation, modulation variation, or both. Furthermore, the 
image contrast may be enhanced by measuring the response of a sample to 
a first excitation light without the agent to provide a baseline (the "absence" 
case), and then measuring the response to a second excitation light after 
20 introduction of the agent (the "presence" case). Data corresponding to these 
two responses is compared to evaluate the contrast capability of agent. The 
first excitation light wavelength may be selected to stimulate intrinsic 
fluorescent response of the tissue at the same wavelength expected to 
stimulate agent fluorescence. Alternatively, the first excitation light 
25 wavelength may be the same as for the fluorescent light emitted by the agent 
to enhance separation of intrinsic tissue fluorescence from fluorescence of 
the agent. Also, multiple comparisons may be performed using different 
wavelengths to better evaluate the influence of the contrast agent. 

In another embodiment of the present invention, the photon fluence 
30 equation and Jacobian estimation process is adapted to determine a map of a 
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designated fluorophore uptake concentration. For this embodiment, a first 
map of chromophore adsorption coefficients Pax— »c scattering 
coefficients m's are determined in the absence of the designated fluorophore 
by estimating the chromophore adsorption coefficient Max-+c ^""^ scattering 
5 coefficient p'g at each grid point j in place of the yield and lifetime estimates. 
Diffusion equation (1) for <Dx(r, cj) may be used in conjunction with modified 
Jacobian equations (7) and (8) to create this first map. The modification 
substitutes the chromophore adsorption and scattering coefficients in place of 
the yield and after adaptation to accommodate these new characteristics as 
10 follows: 
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The elements of the four J acobian matrices employed. JCMx. Ma ). HM^ ^ ), J(Q^ ,x , 

respectively. Updates to the absorption and scattering map were conducted to minimize the 
merit function x - 

n3i,^inai==i^^ Cm J I J ^^^^ 

where n, = Sk and = Di, 

After generating the first map, the designated fluorescent contrast 
agent is introduced, and the total adsorption coefficient determined by 

substituting Max '"^ P'^^^ ^'ax-^c equations (9H11) to obtain a second 
map of the total adsorption coefficient. Noting that Max" l^ax-^m"^ '-'ax-^c* 
and that the uptake of the fluorescing contrast agent is directly proportional 
to Max— >m' uptsl^e concentration may be mapped by determining a difference 
between the adsorption coefficient variations for the first and second maps. 
This "difference map" may then be used to generate an image corresponding 
to the uptake concentration. 

Another alternative embodiment measures the emission responsive to 
each of a number of light source modulation frequencies f. The total number 
of different frequencies employed is designated Mf. To obtain this additional 
data, an iteration of loop 220 is performed for each frequency f indexed to m. 
The number of sources, Sk and detection sites Di are indexed to k and i, 
respectively. This additional data may be used to enhance imaging results 
obtained with system 1 10 or to permit reduction of the number of detection 
sites or excitation source sites in the evaluation. A representative merit 
function corresponding to this additional data is given in equation (12) as 
follows: 



Mf Sk Di 

(1/MO Z(1/Sk) Z(1/Di) Z[((Mobs)i-(MrTn)i)/aM]2^[((eobs)i-(9m)iV^e]2(12) 
m=1 k=1 1=1 
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Besides fluorescence yield and lifetime, the multi-frequency method can be 
employed to map other optical characteristics of interest. Besides a 
sinusoidally modulated light source, the present invention may be adapted to 
operate with a pulsed or other time-varying excitation light source in 
5 alternative embodiments. 

FIG. 15 depicts an optical system 410 of another embodiment of the 
present invention. This system includes modulated light source 420 with 
laser driver 422, operatively coupled laser diode 424, and reference 
frequency generator 426. Source 420 is configured to deliver modulated light 
10 to tissue phantom 400, and the emitted light from the phantom is focused 
onto a gain modulated image intensifier with 50mm lens 432 through filter 
433. Filter 433 may be a bandpass or low pass arrangement selected to 
isolate to a selected wavelength. Typically filter 433 is configured to pass at 
least the anticipated fluorescent emission wavelength and may additionally or 
15 alternatively be configured to pass the excitation light wavelength. Intensifier 
430 includes a photocathode face, which converts photons to electrons, a 
Multi-Channel Plate (MCP) which multiplies the electronic signal by avalanche 
multiplication, and a phosphorescent screen, which converts electrons into an 
optical image. Preferably, intensifier 430 is a fast intensifier, of the variety 
20 manufactured by Litton Electronics. Inc.. which enables modulation by 
applying a DC bias and an RF signal from amplifier 428 between the 
photocathode and the MCP. For this example, the modulation of the image 
from intensifier 430 is phase-locked to the laser diode 424 by a 10 MHz 
output signal from synthesizer 426. By modulating the laser diode 424 and 
25 the image Intensifier 430 at the same frequency, a steady-state image results 
on the phosphor screen. U.S. Patent No. 5,213.105 to Gratton et al. provides 
additional background concerning certain aspects of this technique. The 
image from the phosphor screen is focused through interference filter 433 on 
a Charge Coupled Device (CCD) camera 434 via 150-mm macro lens 436. 
30 Camera 434 has a 512x512 array of CCD detectors configured to provide a 
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corresponding pixelated image. Camera 434 is operatively coupled to 
processor 460 of a similar configuration to processor 160 previously 
described. 

Following each acquired image, a phase delay between the image 
5 intensifier 430 and the laser diode 424 is induced by stepping the phase of 
the image intensifier 430 to values between 0 and 360 degrees with the 
frequency synthesizer 452 under the control of processor 460. Since the gain 
modulation of image intensifier 430 and laser diode 424 occurs at the same 
frequency, homodyning results in a steady phosphorescent image on 
10 intensifier 430, which is dependent upon phase. Preferably, control between 
synthesizer 452 and processor 460 is obtained by a conventional GPIB 
interface. Images from the phosphorescent screen of the image intensifier 
430 are then gathered at each phase delay. The incremental phase delayed 
images are then used to generate a map of phase-shift and intensity 

15 modulation ratio between the excitation and emitted light from phantom 400. 
By applying interference or appropriate optical filters, the emission light may 
be selectively separated from the excitation light and measured. Camera 434 
output may be processed by processor 460 using process 210. 

In other embodiments, a wide area illumination source is preferred to 

20 provide a larger, more uniform front illumination in a reflective geometry. This 
illumination approach facilitates faster imaging of multiple sights and a more 
natural physical correlation between photon migration images and pathology. 
Also, a camera which has a tapered fiber optic coupler from the image 
intensifier to the CCD array is envisioned to increase the efficiency of light 

25 coupling from the intensifier to the CCD array and reduce the physical size 
and weight of the imager. 

The present invention will be further described with reference to the 
following specific Examples 1-8. It will be understood that these examples 
are illustrative and not restrictive in nature. Examples 1-4 involve the 

30 computer simulation of the process 210. Simulations of this kind, including 
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the simulation of tissue, are an acceptable means of demonstrating 
fluorescent spectroscopic imaging perfonmance to those skilled in the art. 
Examples 1-3 use simulated values obtained by solving the diffusion 
eauations (1) and (2 ) for 8^ and under the conditions of table 2 as 
.fpUows: 



Case 






or 


T 




Gatis.<rian 


Gaussian 










(background) 


(background) 


Noise 
in log of AC 


Noise 
in phase 




(mm-') 


(mm-l) 


(iiun-^) 


(ns) 


(mm-^ 


intensity 
<tm 


ce (degrees) 


5.1 


0.0 


0.0 


1.0 


10.0 


1.0 X 10-5 


0.01 


0.1 


5.2 


1.0 X 10-3 


0.0 


1.0 


1 0.0 


1.0 X 10-5 


0.01 


0.1 


5.3 


0.0 


0.0 


1.0 


10.0 


1.0 X 10-5 


0.01 


1.0 



m Table 2 

I The examples simulate tissue phantom 300 of FIG. 3 having a 100mm 

I diameter Values of Q^r^ and were computed at each of the D1-D20 
m detection sites of FIG. 3 in response to the 4 modulated light sources S1-S4 
located at the periphery. The excitation light modulation frequency f was 
simulated at 150 MHz. Diffusion equations (1) and (2) were solved to provide 
80 simulated values of 9^ and corresponding to the various 
combinations of detection and source sites (Sk * Di = 4x20 = 80). Gaussian 
25 noise with a standard deviation of 0.1 degrees (or a liberal 1 degree) in 9^ 
and 1% in were superimposed on the diffusion equation solutions. 
Adapted MUDPACK routines were used to solve the diffusion equations (1) 
and (2) on a SunSpardO computer. These obtained data sets were used as 
simulated input data to process 210 for examples 1-3. The results are shown 
30 in tables 3 and 4 are as follows: 
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Case 


Area, object 1 
(mm^) 


Location, object 1 
(x,y), (mm, mm) 


Area, object 2 
(nmi2) 


Location, object 2 
(x,y).- (mm, mm) 


5.1 


706.0 (expected) 
742.2 (obtained) 


(60,60) (expected) 
(60.8,58.5) (obtained) 


not applicable 


not applicable 


5.2 


706.0 (expected) 

703.1 (obtained) 


(60,60) (expected) 
(59,4,58:3) (obtained) 


not applicable 


not applicable 


5.3 


314.1 (expected) 
381.0 (obtained) 


(32.3,67.7) (expected) 
(34.0,67,7) (obtained) 


314.1 (expected) 
342.0 (obtained) 


(67.7,32.3) (expected) 
(65.0,35.0) (obtained) 



Table 3 



Case 


"^V^^x^m (object) 
(mm*^) 


X (object) 


5.1 


1.0 X 10-3 (expected) 
0.93 X 10-3 (obtained) 


1.0 (expected) 
L03 (obtained) 


5,2 


1.0 X 10-3 (expected) 
0.8 X 10-3 (obtained) 


1,0 (expected) 
0.7 (obtained) 


5,3 


(top left object): 
1.0 X 10-3 (expected) 
2 X 10-3 (obtained) 
(bottom right object): 
2,0 X 10-3 (expected) 
1.8 X 10-3 (obtained) 


(top left object): 

1.0 (expected) 

4.1 (obtained) 
(bottom right object): 

2.0 (expected) 
3.5 (obtained') 



Table 4 




Example 1 

Example 1 reconstructs fluorescent yield and lifetime with no 
absorption due to non-fluorescing chromophores. To simulate the 
5 experimental data for this example, the fluorescent yield, (nPax-^m^j' ^® 
background and the heterogeneity 302 were chosen as 1 x 10"^ mm"'' and 1 
x 10-3 mm-T respectively and the fluorescence lifetime. (T)j, for the 
background and the heterogeneity 302 chosen as 10 ns and 1 ns 
respectively. During the execution of loop 220. no a priori knowledge of either 
10 the heterogeneity 302 location or the background fluorescence properties was 
assumed and a uniform guess of 1x lO'^ mm*'' and 10 ns was given for the 
fluorescence yield. (nMax-^m^j' lifetime. (T)j. respectively. Convergence 
was achieved in less than 50 iterations of Loop 220 (computational time on a 
SunSpardO: 2 hours) for a two dimensional 17x17 grid. The average 
15 values of HMax— ^m ^ 9''*^ points which occupy the simulated object 

converge within 50 iterations to HMax^m " ^-^^ ^ ^^'^ mm"^ and x =1.03 ns 
are illustrated in FIGS. 8 and 9. respectively. FIGS. 10 and 11 illustrate the 
reconstructed images from the mapped values of ri|Jax-.,.rn ['t^'^^"^] ^nd x [ns], 
respectively, and are representative of the expected images. The images 
20 were smoothed by interpolation in examples 1-3 to remove spurious points 
which had unphysically high values, but were surround by values within a 
physically achievable range. These spurious values were replaced by the 
average background fluorescence yield and lifetime obtained from simulation 
of loop 220. 

25 The average values of qPax— m S*"'^ points which occupy the 

simulated background converge within 50 iterations to 9 x lO'^ mm"''. The 
value of the background converges to 5.4 ns. The dependence of the final 
images on the choice of the initial guess was examined by providing an initial 
uniform guess of 1 x 10"^ mm*^ and 10ns for (HMax— ^m^j' lifetime, (x)j, 

30 respectively. This resulted in similar images to those obtained in FIGS. 10 




and 11 . 

The location of heterogeneity 302 was identified as consisting of ail the 
grid points with nPax-^m ^^'O^er than 35% (arbitrarily chosen) of the peak 
value of the HMax-^m '^^^ average of the coordinates of all the 

5 identified object grid points was the position (60.8, 58.5) which is close to 
position (60, 60) that was used to simulate the experimental data. As listed in 
Table 3, the area of the heterogeneity based upon our arbitrary definition for 
identification was 72 mm^, close to that used to generate our simulated 
experimental data. 

10 

Example 2 



Example 2 reconstructs fluorescent yield and lifetime with a simulated 
chromophore absorption configured to mimic tissue. The same hidden 
15 heterogeneity as well as optical parameters and simulation equipment were 
used as described in Example 1 except that a uniform background 
chromophore absorption coefficient, Max— ^ xlO"^ mm"^ was used to 
generate the simulated experimental data. While excitation light propagation 
was not employed for image reconstruction, we considered this optical 

20 property known to estimate the best possible performance for inverse image 
reconstruction under physiological conditions. The two-dimensional 
reconstructed spatial map of the fluorescence yield, (RMax— ►m^j ["^"^"^J* 
lifetime. (x)j [ns], are shown in FIGS. 12 and 13, respectively. As shown in 
Table 3, the mean value of location of the object according to our criterion 

25 based on nPax-H^m occurred as position (59.4, 58.3) consistent with the 
conditions used to simulate the experimental data. The dimension of the 
heterogeneity based upon our arbitrary definition for identification (all grid 
points with HPax— ^m higher than 35% of the maximum) were 703 mm^ which 
is close to that used to generate our simulated experimental data. The 

30 average values of nPax-^m '"^ ^^^^ points which occupy the 
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simulated object converge within 50 iterations to the values of r]Pax— i-m "0*8 
x10"3 mm"'' and x = 0.7 ns consistent with the values used to generate the 
simulated experimental data (see Table 3). The average values of HMax— ►m 
and T in the grid points which occupy the simulated background converge 
5 within 50 iterations to values similar to that reported for Example 1. 

Example 3 

Example 3 ' nulated two hidden heterogeneities in the tissue phantom 
(not shown in FIG. 3). In this case, the same optical parameters were used 
10 as described in example 1 except that the fluorescence yield HMax— •►m 

objects 1 and 2 was chosen as 1 x lO'^mm"'' and 2 x lO'^mm"'' respectively 
and lifetime t for the heterogeneities chosen as 1 ns and 2 ns, respectively. 
A 33 X 33 grid was employed instead of a 17 x 17 grid. An image 
corresponding to the mapping of yield is depicted in FIG. 14. 

15 

Example 4 

Example 4 demonstrates the unexpected advantage of utilizing a 
fluorescent contrast agent with a fluorescence lifetime within an order of 
magnitude of the mean time-of-flight of the interrogating photons. This 

20 example compares by computational simulation, the contrast offered by a 
phosphorescent agent with a lifetime of about 1 millisecond to a fluorescent 
agent with a lifetime of about 1 nanosecond. Referring to Fig. 16. this 
simulation assumes a circular tissue phantom 500 with an embedded 
heterogeneity 502. A single source SI is indicated and detectors D1 through 

25 D1 1 are spaced about halfway around the circular periphery at generally 
equal intervals. For this comparison ujt = 1 for both the heterogeneity its 
surroundings. An uptake of the heterogeneity of 100 times the surroundings 
(background) was assumed with a corresponding absorption coefficient of 

fluorophores in the heterogeneity set to 0.1 cm*^ versus 0.001 cm"^ for the 
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surroundings. The absorption coefficient for nonfluorescing chromophores 
was set to 0.001 cm""" for both the heterogeneity and surroundings. 
The graph of Figs. 17A plots measurements of phase shift (vertical axis) 
versus angular detector location about the circumference of the tissue 
5 phantom 500 (horizontal axis). The lines with open symbols show phase shift 
variation of a contrast agent with a lifetime of 1 ns with the detector position 
and changes in location of a heterogeneity containing the contrast agent 
(different open symbol shapes). The closed symbols correspond to phase 
shift versus detector position for a contrast agent in-the heterogeneity having 
10 a lifetime of 1 ms. The graph of Fig. 17B compares these contrast agents in 
terms of Modulation (vertical axis) versus angular detector location (horizontal 
axis) and location of the agent-containing heterogeneity (different open 
symbol shapes). These illustrations indicate that as the lifetime increases, 
sensitivity of the contrast agent to spatial differences is reduced, and further 
15 point to the development of fluorescent contrast agents with intrinsic 

relaxation, kinetics (as characterized by fluorescence lifetime) within an order 
of magnitude of the photon migration times (i.e. times-of-fiight) for imaging 
based on the behavior of multiply scattered light. 

Figs. 18A and 18B are calculated simulations using the same phantom 
20 tissue and parameters, comparing lifetimes of 1 ns and 1 ms, respectively, as 
a function of position of the detectors and location of a heterogeneity 
selective to the corresponding contrast agent. The vertical axis of Figs. 18A 
and 18B represents phase contrast AG. Phase contrast AQ is the difference 
between the phase shift in the presence of the heterogeneity and the phase 
25 shift in the absence of the heterogeneity, (Ae = Gpresence - 9absence)- 

Figs. 18C and 18D are calculated simulations using the same phantom tissue 
and parameters as Figs. 18A and 188 to compare the same contrast agent 
lifetimes in terms of modulation contrast, AM (vertical axis). Modulation 
contrast AM is provided as the ratio of the modulation ratio (AC/DC) of the 
30 detected light in the presence of the heterogeneity to the modulation ratio 
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(AC/DC) of the detected light in the absence of the heterogeneity, ( AM = 
Mpresence^Mabsence)- ^^r all of Figs. 18A-18D. the horizontal axis 
corresponds to the detector number and the different line styles correspond to 
different positions of the contrast agent-bearing heterogeneity. 

5 

Example 5 

The conclusions of the simulation of Example 4 have further been 
empirically demonstrated by the experimentation of Example 5. The 
experimental equipment set-up for Example 5 is comparable to system 110. 

10 A tissue phantom is prepared by filling a cylindrical Piexiglas container having 
a 20 cm diameter and a 30.5 cm height with a Q.5% Intralipid solution 
(supplied by Kabi Parmacia. Clayton, NC). A heterogeneity is provided by 
placing a cylindrical glass container with a 9 mm inner diameter in the 
Piexiglas container and filling the glass cylinder with the intralipid solution and 

15 a contrast agent. The position of the heterogeneity within the Piexiglas 

container is adjusted with an x-y translation stage model number PMC200-P 
supplied by Newport of Irvine. CA. 

Example 5 experimentally confirms the phase contrast simulated in 

Example 4 by comparing phase contrast A9 with a Ru(bpy)3^'^ 

20 phosphorescent contrast agent in the inner glass container (Fig. 19B) to a 
fluorescent contrast agent. Indocyanine Green (ICG) (supplied by ACROS 
Organics, Fairlawn. NJ). in the inner glass container (Fig. 19A). The contrast 
agents where added to the intralipid solution in the inner glass container (the 

"heterogeneity") to simulate a 100:1 uptake ratio. Ru(bpy)3^* has a lifetime 

25 on the order of microseconds. ICG has a lifetime of about 0.58 ns. Notably, 
by comparing Figs. 19A and 19B, the phase contrast (vertical axis) provided 
by a fluorescent agent with its smaller lifetime, is substantially greater than the 

phase contrast provided by the longer-lived phosphorescent agent Ru(bpy)3^* 
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as a function of detector number (horizontal axis) and heterogeneity location 
(different line styles). 

Example 6 

5 The experimental equipment set-up for Example 6 is comparable to 

Example 5, except a single source and a single detection point were utilized. 
The source and detector were placed along the circumference a few degrees 
apart and the inner container was generally positioned along the midline 
defined between the source and detector. The x-y stage was used to adjust 
10 the position of the inner container along this midline to observe corresponding 
changes in phase shift 9 and amplitude. 

In Example 6, the response of two different fluorescent contrast agents 
ICG and 3-3'-Diethylthiatricarbocyanine Iodide (designated "DTTCT herein 
and supplied by ACROS Organics. Fairlawn, NJ) to an intensity-modulated 
15 excitation light having a wavelength of about 780 nm was detected. The 
excitation light was modulated at 80 MHz and 160 MHz in different trials 
corresponding to lines with different symbol shapes. ICG is an agent 
approved for hepatic and retinal diagnostic testing with a measured lifetime of 
about 0.58 nanoseconds and DTTCI is a common laser dye with lifetime of 
20 about 1.18 nanosecond. 

The phase shift and modulation ratio of the tissue phantom in the 
absence of the heterogeneity was measured to provide the "absence" case 
needed to calculate phase contrast A9 and modulation contrast AM. Next, an 
ICG contrast agent was prepared by adding about a 2.0 fjmole ICG 
25 concentration to a 0.5% intralipid solution in the inner container. The ICG 
sample was then exposed to excitation light from the source and the 
response detected. This detection included measurement of absorption at a 
wavelength of 780 nm and fluorescence at 830 nm. The resulting phase 
contrast A9 at the absorption wavelength (open symbols) and for 
30 fluorescence wavelength (closed symbols) for the ICG sample was plotted on 
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the vertical axis of the graph provided in Fig. 20A with the horizontal axis 
showing the relative position of the heterogeneity ("object") in centimeters as 
it is moved toward the detector and source along the midline. The resulting 
modulation contrast AM at the absorption wavelength (open symbols) and at 
5 fluorescence wavelength (closed symbols) for the ICG sample was plotted on 
the vertical axis of the graph provided by Fig. 20B with the horizontal axis 
showing the relative position of the object in centimeters as it is moved toward 
the detector and source along the midline. 

After ICG sample was tested, a 4.2 Mmole concentration of the DTTCI 
10 contrast agent was added to the 0.5% intralipid solution in the inner container 
to provide a DTTCI sample. The different concentration of the ICG and 
DTTCI contrast agents were selected to provide a fluorescent cross-section 
that is generally the same for both the ICG and DTTCI samples. The 
resulting phase contrast AS at the absorption wavelength (open symbols) and 
15 at the fluorescence wavelength (closed symbols) for the DTTCI sample is 
plotted on the vertical axis of the graph provided in Fig. 20C with the 
horizontal axis showing the relative position of the object in centimeters as it 
is moved toward the detector and source along the midline. The resulting 
modulation contrast AM. at the absorption wavelength (open symbols) and for 
20 fluorescence wavelength (closed symbols) for the DTTCI sample is plotted on 
the vertical axis of the graph provided in Fig. 20D with the horizontal axis 
showing the relative position of the object in centimeters as it is moved toward 
the detector and source along the midline. 

For both samples, the fluorescence decay process is single 
25 exponential, showing one lifetime, but the analysis and approach can be 
extended to dyes and contrast agents with more than one lifetime. Upon 
comparing the fluorescent phase and amplitude modulation generated by the 
two fluorescent contrast agents, the impact of fluorescence lifetime x over 
absorption may readily be observed. Indeed, it has been found that 
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substantial contrast is present when the uptake is only 10:1 over the 
surroundings or background. 

Example 7 

5 For Example 7, a phantom tissue is prepared by placing a tissue- 

mimicking Intralipid solution in a Plexiglas container. The excitation light 
source transilluminates the tissue phantom from the rear along a straight-line 
distance of about 8 centimeters. An image intensifier/CCD detection 
arrangement was utilized to detect the response. The experimental set-up for 

10 Example 7 was comparable to system 410 illustrated in Fig. 15. 

Embedded within the middle of the Plexiglas container tissue phantom 
were two micromolar intralipid solutions of 0.5 ml in separate containers each 
having a different fluorescent contrast agent. One vessel included an ICG 
contrast agent and the other vessel included a DTTCI contrast agent. 

15 Measurements of the fluorescent phase-shift, AC amplitude, DC intensity, 
and modulation (AC/DC) were conducted across the front of the phantom 
tissue in response to a 100 MHz modulated excitation light at 780 nm. Figs. 
21A-21D are two-dimensional images depicting the spatial variation of the 
prtase-shift. AC amplitude, DC intensity, and modulation measurements; 

20 respectively, in terms of a corresponding gray scale. These images confirm 
the variation in contrast with differences in fluorescence lifetime and the 
parameter being measured. 

Example 8 

25 Example 8 is a live tissue study of mammary tissue from a dog, Sugar 

Limburg, which was a miniature poodle (age 10 years and weight 12.5 lbs.). 
An in vivo image of the right fifth mammary glad was taken after an in vivo 
injection with 1.3 cc of a 5% concentration of ICG fluorescent contrast agent. 
Interrogation was performed with an experimental set-up comparable to 
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Example 7, with an excitation light wavelength of 789 nm and detection at a 
830 nm wavelength. The modulation frequency was 100 MHz. 

A frozen section of the right fifth mammary revealed two dark spots 
approximately 1 cm deep from the tissue surface which were histologically 
5 classified as reactive regional inguinal lymph nodes with no evidence of 

metastatic spread. The remaining tissue was classified as lobular hyperplasia 
with no evidence of tumor. Figs. 22A-22D are 2-dimensional images 
depicting spatial variation in terms of in vivo measurement from the emission 
light for modulation phase, modulation ratio, average intensity, and 
10 modulation amplitude, respectively, relative to corresponding gray scales. 
The white line in the image of Figs. 22A-22D is the air-tissue boundary. The 
phase shift in the tissue is small due to the small amount of mammary tissue 
present. Also, the modulation ratio contains a large amount of background 
noise. However, the two light spots at the bottom of both the average 
15 intensity and the modulation amplitude are believed to correspond to an 
increased uptake of ICG inside the enlarged lymphatic tissue. 

All publications and patent applications cited in this specification are 
herein incorporated by reference as if each individual publication or patent 
application were specifically and individually indicated to be incorporated by 
20 reference; including United States Patent Application Serial Numbers 
60/039,318 filed 7 February 1997 and 08/702,060 filed 23 August 1996. 
While the invention has been illustrated and described in detail in the 
drawings and foregoing description, the same is to be considered as 
illustrative and not restrictive in character, it being understood that only the 
25 preferred embodiments have been shown and that all modifications that come 
within the spirit of the invention are desired to be protected. 



